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Abstract 



We study the effects of gluon interference in the production and semi-leptonic 
decay of a tt pair above threshold at the Next Linear Collider (NLC). We calcu- 
late all matrix elements to next-to- leading order and use the resulting expressions 
for the development of a Monte Carlo event generator. Our results show effects 
at the level of 10% in differential cross-sections. We thus extend previous results 
obtained by analytical means in the soft-gluon limit. 
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I. INTRODUCTION 



The discovery of the top quark |],|2| has opened up a new avenue for the exploration of new 



physics. The large mass of the top {m t = 176 ± 8 ± 10 GeV (CDF) 0; m t = 199+^ ± 22 
GeV (DO) pi) and, consequently, its large decay width (r t = 1.57 GeV 0), permit the use 
of perturbation theory, because the top quark is expected to decay before it has a chance to 
hadronize. It is unlike the other (lighter) quarks, which form resonances that are outside the 
realm of perturbation theory. It is therefore possible to perform a precision analysis of top 
quark production and decay, which facilitates the identification of effects attributable to new 
physics beyond the Standard Model ||. 

The best laboratory for such explorations is an e + e~ collider, because the initial state is free 
of color effects ||. Moreover, it is advantageous to have polarized beams for a more efficient 
study of the electroweak properties of the top quark. Such a machine (the Next Linear Collider 
(NLC)) will hopefully be built in the not-so-distant future. Detailed studies of processes through 
which one may extract useful information have already been performed 0. Gluon radiation 
contributes at the 10% level, in general. Although it does not modify inclusive cross-sections, it 
has a non-negligible effect on various final-state distributions. Therefore, it has to be included 
in a precision study of top quark physics. In studying gluon radiation, one cannot ignore the 
effects of gluon interference. The latter depends upon the value of T t and may be an important 
tool in the study of this parameter, in addition to precision studies [0]. 

The decay width of the top quark is saturated by the decay mode t — > W + b, because in the 
Standard Model \V t b\ ~ 1. The W boson, in turn, may decay into quarks or leptons. Therefore, 
there are four possible decay modes in tt production and decay. The branching ratio for the 
channel in which both W bosons decay into leptons is about 10%. This is the channel we study 
here, because it gives the cleanest signal. We shall calculate all amplitudes analytically, with 
the aid of the symbolic manipulation package FORM M. The resulting expressions will then 
be used for the development of a Monte Carlo event generator. The generator is based on the 
C++ program written by Schmidt for the study of the same process M, 

e + e~ -> 7* , Z* -> tt -> W + b W~b -> tvb t~vb (1) 

We have added gluon interference effects to the generator. The results from our simulations 
are in agreement with analytical results obtained by Khoze, et al. |TD| in the soft gluon regime. 
They are also similar to the results obtained by Peter and Sumino [|ll|] in the threshold region, 



where color Coulomb effects dominate IT2 



Our paper is organized as follows. In Section |J, we outline the tree-level analysis in order to 
introdce the basic tools and establish the notation. In Section [TTI|, we discuss loop amplitudes at 



o(a s ). In Section [TV], we introduce real gluons and compute the amplitudes for gluon radiation. 
In Section M, we discuss the development of the Monte Carlo event generator based on the 
analyical results obtained in the previous sections. We also present numerical results. Finally, 



in Section NIL we present conclusions and discuss future directions. 
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II. TREE-LEVEL ANALYSIS 



At tree level (Fig. |l|), the amplitude for tt production and decay factorizes. With polarized 
e + e~ beams, one can perform a clean helicity analysis, because the spins of the top and anti-top 
quarks get transferred to the final products. It is then advantageous to express the amplitudes 
in terms of helicity states ||. To this end, one defines form factors for the production and 
decay stages. 

The form factor for tt production is the neutral-current coupling (7-t-t and Z-t-t vertices) 

= erf(F% + F}frf) + l^f V K{F™ + F 7 2 l 7 5 ) , (2) 

where in the Standard Model the only non- vanishing coupling constants at tree level are 

F] V = Qt = - , Fz v = - - tan6 w , F\ A = - . (3) 

' 3 2sm{20w) 3 2sm(2tV) 

The magnetic dipole coupling (F 2V ) receives loop corrections at o(a s ). The electric dipole 
coupling (F 2A ) is CP violating and remains zero to o(a 2 s ). 

The t — > W + h decay vertex is described by the form factor (top-charged-current coupling) 

M» w = ^(FwPl + FwPr) + ^^f Vk ^ F % P L + F wPr) , (4) 

where Pr,l — =t 7 s ), and similarly for the i — > W~b decay vertex. In the Standard Model 
at tree level, F^ = 1 and all other coupling constants vanish. 

Finally, the W + — > t + v decay is described by a left-handed current coupling of strength 
«S and similarly for W~ -> tv. 

Scattering amplitudes can be conveniently expressed in terms of helicity angles 0. There 
is a set of helicity angles for each vertex. Different sets are defined with respect to different 
Lorentz frames. This assumes the factorizability of the amplitudes. Since we want to study 
gluon interference, we have to consider non-factorizable amplitudes at one-loop level. Therefore, 
we will use a single Lorentz frame (laboratory frame). To simplify the expressions, we introduce 
polarization vectors for the W ± bosons defined with respect to the b (b) quarks (we set = 
in all matrix elements), respectively, 



Pw+ ■ Pb 



(w $ y = -^^--™-rf ( 5) 

v 7 m w p w + - pb 

and similarly for (W~) fl (only one transverse polarization of each W contributes to matrix 
elements). Longitudinally polarized W bosons are very important because the study of longi- 
tudinal modes probes the electroweak symmetry-breaking mechanism. 
The tree-level amplitudes can be written 
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T 7>Z (L, W-, W + ) = (F^ z F^ z}a T^(L, W~ , W + ) - m 2 t F^ z F^ z}a T^(L, W~, W + )) W + W~ 

(6) 

where we chose left-handed polarization for the electron beam, for definiteness. F^ z and FP z ^ 
are appropriate form factors given in terms of the form factors in Eqs. (0), (^), and (0j), and 

T^(L,W~,W + ) = (b R \ W fc\ei)(el\ AW + *\b L ) 

T^(L,W-,W + ) = (e+| W + *\b L )(b R \ W~\el) (7) 

The factors 

w + = (41 w + Wl) , w- = {v R \ w~yi) (s) 

represent the leptonic decays of the W ± bosons, respectively. The differential cross-section for 
left-handed polarized beams at tree-level is 



da tree (L) = |£ (T 7 (L, W~, W + ) + T Z (L, W~, W + )) (9) 
where we sum over the polarizations of W ± . 

III. VIRTUAL GLUONS 
A. General 

Loop diagrams at o(a s ) include a pentagon diagram (gluon exchange between the b and b 
quarks; see Fig. $(f)), which renders the calculation of o(a s ) effects cumbersome. To systemat- 
ically calculate these effects, observe that all loop diagrams contain one or more internal t or t 
legs, which are rapidly decaying quarks. To take advantage of this fact, we define the quantities 

Vi = (Pi ~ m t)/ m t , V2 = (Pt - m 2 )/m 2 ■ (10) 
Due to the presense of propagators 

1 



i>-m t - iT t /2 ' 

where is the momentum of t (t), the integral over phase space contains factors 

drji f dr] 2 



(11) 



r, 2 + Y 2 /m 2 t J ri + Y 2 t /m 



(12) 



showing that 771 and 772 are small parameters, o(T t /m t ). Physically, this means that t and t are 
produced nearly on shell. 
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The rest of the amplitude is a function of r\\ and r)%. Since r]\ and r\ 2 are small, it is 
advantageous to expand the amplitude in these parameters, 



•A(Vi,V2) = A 



OA 

fj! =7)2=0 (jf]l 



V2=0 



OA 

dr] 2 



(13) 



where A needs to include both virtual and real gluon contributions (soft and collinear), in order 
to avoid singularities, 



A = A 



virtual 



+ A soft + A 



collinear 



(14) 



The first term in Eq. ( |T3"D is the set of triangle diagrams (Fig. ^a-cj). This set of diagrams is 
gauge-invariant if we set r/i = rj 2 = 0. They represent vertex (form-factor) corrections |J. The 
rest of the terms represent gluon interference. To calculate them, notice that if we set r] 2 = 0, 
then we are placing the top quark on shell, thus introducing a cut. The diagrams contributing in 
this case are those containing a gluon which does not cross the cut (Fig. ^( a,b,d) and Fig. Q( a-c)) . 
This set is gauge-invariant at the point rj 2 = 0. Notice that the pentagon diagram (Fig. ^(f)) 
is not included in this set. Their contribution to the amplitude is 



Ai=A 



V2=0 



OA 

vi=V2=o arji 



(15) 



r?2=0 



Similarly, if we set 771 = by cutting the t line, we obtain a set of diagrams which is gauge- 
invariant at the point r/i = (Fig. ^(a,c,e) and Fig. ^(b-d)). Their contribution to the amplitude 
is 



A 2 =A 



Vi=0 



A 



dA 

Tj! =7/2=0 0*772 



+ o(vl 2 ) 



(16) 



The pentagon diagram (Fig. ^(f)) contributes to o(r]f 2 ). 
The o(a s ) amplitude (Eq. (|13D ) can thus be written 

A( V i, r} 2 ) = A(0, 0) + (A - A(0, 0)) + (A 2 - A(0, 0)) + o( V l 

We shall calculate each term separately. 



(17) 



B. The limit of on-shell t and t 



In this case, the amplitude factorizes. There are three triangles (Fig. ^(^a-c^)) contributing 
and have already been calculated 0. Here, we express them in terms of W ± polarization 
vectors (Eq. (|5|)), mainly for completeness. For definiteness, we present results for left-handed 
beam polarization. Right-handed polarization can be obtained by interchanging e + and e~ . 

The o(a s ) correction to the tt production vertex is 



C F (gn 



e\2 



d*-*k C HZ} u(L,W-,W+) 



(27r) 4 - 2e k 2 (k 2 - 2k ■ P t)(k 2 + 2k ■ p t ) 
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where 

C HZ} u(L, W~, W+) = F^ L z F^ z}t tCV(L, W~,W + ) - m 2 t F^ z F^ z}t £ ( ?' '(L, W~, W + ) (19) 

Cf{L, W-, W + ) = 2p t ■ Me L \ i t W + *\b L )(b R \ W~ fc\e+) 
+ (e L \ £ tfaW + *\b L )(b R \W- AK) 
-{el\A t W + *\b L ){b R \W' AtA^V) 
- (41 ¥ A W + *\b L )(b R \ W~ A ¥\el) (20) 

C%\L,W-,W + ) = -2p t - Pl (e+\ W + *\b L )(b R \ W'\e~ L ) 
-(e+| ^- t W + *\b L }(b R \W~\e- L ) 
+ (e + R \W + *\b L )(b R \W- A^l) 
+ (el\ ftW + *\b L )(b R \W- ¥\e + R ) (21) 

These expressions can be further simplified by employing helicity states for the on-shell t and 
t quarks. We shall not do it here, because we are interested in the more general case of off- 
shell quarks. Notice that we obtain the correct infrared limit (after we subtract ultraviolet 
divergences) 

° F{9fi } J (2^ ¥ U 2 - 2k - Pt - ~ P + 2k- ih) T ^(L,W,W) (22) 
The o(a s ) correction to the t — > W~b decay vertex is 

e . 2 r d^k C HZ} t(L,W~,W + ) 



J (27r) 4 - 2e k 2 {k 2 -2k-p- b ){k 2 -2k-pt) y ' 
where 

C HZ}i (L, W~,W + ) = F^ z F^ z}t ^(L, W~, W+) - m 2 t F^ z F^ z]t ff\L, W~,W + ) (24) 

Ci 1] (L,W-,W + ) = (e- L \ ^W + *\b L )(b R \IA\e^ (25) 

Cf(L, W',W+) = (l R p\eZ)(e%\ W + *\b L ) (26) 

J = 2(-pt + k) ■ p- b W ~ 2Pb ■ W~ j6+ i>- t $ W~+ ¥ W~ P (27) 

Again, the correct infrared limit is obtained, similar to Eq. (P2"|). The o(a s ) correction to the 
t — > bW~ decay vertex can be obtained from the t — > bW + decay vertex by conjugation. 
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C. The limit of on-shell t 



Tj), (H) and (H)) 



Here we calculate the virtual gluon contribution to A\ (Eq. (|15|)). The vertex corrections 
we calculated above need to be modified to include the o(rji) effects of an off-shell t quark. We 
obtain additional contributions (Eqs. 



(28) 



ACit\L,W-,W + ) =m 2 Me- L \ ^ t W + *\b L )(b R \ W~{*t+ J0I4> 
-m 2 M4\ ^ t W + *\h)(b R \W~\el) 



ACP(L,W-,W + )=2(-p- b + k)-W-m 2 t r ]l (e- L \ t W + *\b L )(b R \e+) (29) 

and AC^(L, W~, W+) = ACf\L, W~, W + ) = 0. 

In addition to corrections to the triangle diagrams of 0(771), we have to consider the box 
diagram in which a gluon is exchanged between b and the top quark (Fig. @(^). Its contribution 
is formally of o(r]i), 

r , ~ 2 f d^k V^ Z (L,W-,W+) 

FWHl i (2vr) 4 - 2e k 2 (k 2 -2k-p~ b )(k 2 -2k- P t + m 2 r ]l )(k 2 + 2k-p t ) 1 J 

where 

V,, Z (L, W~, W + ) = F^F^gV^L, W~,W + ) - m 2 F^ L z F^ z}t ^ 2 \L, W~,W + ) (31) 
VV>(L, W~, W+) = (el\ AW + *\b L )(bM t + ft){- W~(- A+ ¥)\e + R ) 

- H(- A+ V) W~(- A+ WI4)<&*| ¥ A W + *\b L ) (32) 
V^(L,W-,W + ) = (e+| W + *\b L )(b R \ $ t {- fa+ ft) W~\el) 

+ (41 K~ A+ ¥) W~\el)(b R \ W + *\b L ) (33) 



To perform the integral over the loop momentum in Eq. (|30|) , we can use a standard reduction 
procedure implemented via the symbolic manipulation program FORM ||. The only cumber- 
some terms are those with three loop- momentum insertions in Eq. (|3~2"D, 



(e L \AW + \b L )(b R \ ftftW Mr) -(el| ftW~ 04><&h| ft A W + \b L ) (34) 

The first term in Eq. (|3~4j) is readily reduced to a triangle diagram, because ft ft = k 2 cancels 
one of the propagators in Eq. fl3"0p. To reduce the second term, we note that 



¥= 2k -W~ f£- k 2 W~ (35) 
For transverse polarization, we have (Eq. (||)) 
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k - W ~ = T7^ &l M ^ ( 36 ) 

W2p w - ■ p b 

When multiplied by (bn\ ft . . ., we obtain an expression of the form 

. . . ft f b ft ... = ... (2k ■ p b ft - k 2 f b ) ... = ... (- ft(k 2 - 2k ■ p- b ) + (- f b + ft)k 2 ) . . . (37) 
which only contributes to triangles. For longitudinal W~ polarization (Eq. @), 

k ■ W~ = — !— (k ■ pt - Vw ~' Vl k ■ p- b ) (38) 
m w \ Pw- ■ Pt J 

which leads to triangles and a box with one less loop momentum insertions. This is the advan- 
tage of introducing polarization vectors for W . 

D. The limit of on-shell t 



The amplitudes contributing to A% (Eq. flTEp), which include 0(772) effects, can be obtain 
from the expressions we derived above for A\ by conjugation. 

E. Both t and t off shell 

In the case where both t and t are off shell, we need to introduce corrections to the diagrams 
we considered above, as well as the pentagon diagram with a gluon exchange between the b and 
b quarks (Fig. ^/j). These additional contributions to the amplitude are of o(r/^ 2 ) an d will not 
be considered here. We shall only point out that they do not present additional computational 
difficulties. They can be reduced to scalar amplitudes in the same manner. 

IV. REAL GLUONS 

The calculation of amplitudes for real gluon emission (Fig. |3|) follows the same lines as the 
calculation of the virtual contributions. We introduce polarization vectors for the gluon defined 
with respect to the b or b quark, accordingly (cf. Eq.(^)), 

el = 4vi Pb ■ k {hl [r ' ' eR = el (39) 

where k^ is the momentum of the gluon, and similarly for the b quark. In the limit of on-shell 
t and t quarks and the narrow-width approximation, one can clearly distinguish between three 
channels of gluon production (tt production, t decay and t decay, respectively). These channels 
get blurred by gluon interference effects. However, it is advantageous to keep them distinct from 
a programmatic viewpoint. To this end, we split the phase space into three distinct regions. 
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A. The t decay channel 

In this channel we define p^ = p 1 ^- + p^ + k^ and p% = p^ v+ + p%. 77^2 are still given by 
Eq. flTO|). We set r]2 = in the matrix elements, so that the integral over soft and collinear 
gluons contributes to A\ (Eq. [15|). The region of phase space corresponding to this channel is 
defined by 

\r]im 2 t \ < \{pt~ k) 2 - m 2 t \ (40) 

There are three contributions to this channel; gluon emission from b, t and t (Fig. Bl^a-cj). 
They are, respectively, 

R- h = -J—(&I&7&-w*lZF3 F 7& 



-2k -p- b 



7 ■ytt b t 7 "fit 



= (b R \ # w~ ^I4)(ell A W + *\b L ) 

*f = (e+\ W + *\h)(b R \ W~\el) (41) 

* = -2k- p)^ ^ ~ 
t4 ] = (b R \ W~{(- M JO iit-m\ A\e$(el\ AW + *\b L ) 

= (e+| W + *\b L )(b R \ W-{-2e-p;+ ft <fi\el) (42) 

Rt ~ -2k.p- b (-2k.p- t + m 2 r ll ) ^ F ^ " ^ F ^ ) 

= (b R \ W'(- ti+ fl\ei)(el\{(A+ ¥) tit +™ 2 t fi W + *\b L ) 

TZf = (et\{2e-p t + jl W + *\b L )(b R \ VHO (43) 

This set is gauge-invariant in the limit r] 2 = (if we set e M = k^, we obtain 7£& + TZt + TZ t = 0). 
The soft and collinear limits are defined by k° < u and kj_/k° < xq, respectively, where Uq 
and Xq are arbitrary small parameters. Let Cl be the corresponding region of the gluon phase 
space. The contribution to A\ is 
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B. The t decay channel 



This channel can be obtained from the t decay channel by conjugation. Specifically, we 
define p% = p 1 ^- +Pr and p^ = p^+ +Pb + ^ with 771,2 given by Eq. ([11]). We set 771 = in the 
matrix elements, so that the integral over soft and collinear gluons contributes to A2 (Eq. fll6l))- 
We define the gluon polarization vector with respect to the b quark (c/. Eq. ( j39|) ). The phase 
space is defined by \r]2m 2 \ < \(p t — k) 2 — m||. The amplitudes that contribute correspond to 
gluon emission from the t, t and b quarks (Fig. H(7>-(Zj). This set is gauge invariant in the limit 
771 = 0. One obtains expressions similar to Eqs. (|41|) , (|42|) and (|43|) . The amplitude for real 
gluon emission from b is 

k£ ] = (b R \ W~ fc\ei)(el\ A W + \A+ ¥) Ab L ) 

^ = (41 W + *(A+ ¥) Ah)(b R \ W-\e L ) (45) 
The amplitude for real gluon emission from t is 

2k ■ p t + mf i] 2 v ' ' ' ' 

= (b R \ W~ A\ei)(e L \y t ¥)+™ 2 t /} W + *\b L ) 

Tlf ] = (e+\{2e-p t - M)} W + *\h}(b R \ W'Vl) (46) 
The amplitude for real gluon emission from t is 

2k ■ p b (2k ■ p t + m 2 r] 2 ) 



D__ m t^2 / L L (1) 2 L R ( 2 )\ 



nf ] = (l R \ yr-{* AM ¥) - ™ 2 t A\e + R )(eM t - ¥) W + *K) 

Tlf = (e+l W + *\b L )(b R \ W~{2e- P t+ (47) 
contributing to A2 in the soft (A; < u ) and collinear (k±/k° < x ) limits 

C ^ e f I ^3- 2 ! /0 El^ + ^- + ^| 2 (48) 



In (2vr) 3 - 2e A; 

where Q is the region of gluon phase space in these limits. 
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C. The tt production channel 



This channel contains contributions from gluon emission from all quark legs. Since we 
are only keeping terms that are linear in 77^2, we may divide the channel into two parts, 
corresponding to the two channels we considered above. The 0(771) (0(772)) terms are given 
by the 0(771) (0(772)) terms in the t (t) decay channel. The only difference is that in this channel 
we define p% = p 1 ^- + Vi an d Pt = Py/+ + Pb ■ The phase space is restricted by 

1 771m j I < I (pi + k) 2 - m 2 I , \r\ivr? t I < I (p t + k) 2 - m 2 | (49) 



In the limit 771 = 772 = 0, we obtain identical expressions from Eqs. (fll| - f43| ) and Eqs. (45 - 47). 
In the soft limit (k° < ujq), we obtain 

- c ^ 2 U te - w^ T ^ w '- w+) (50) 

which is the correct infrared limit. 



V. THE MONTE CARLO EVENT GENERATOR 

We used the formulas obtained above for the amplitudes to construct a Monte Carlo event 
generator. The generator is based on a program written by Schmidt in C++. Schmidt's program 
generates events at tree level with 100 % efficiency This is achieved by making use of helicity 
states, not only for the final states, but also for the intermediate t and t quark states. The 
amplitude is then expressed in terms of helicity angles which are defined in different Lorentz 
frames for the three vertices, respectively. This construction is carried over to one-loop level, if 
one neglects gluon interference effects, because the amplitude is still factorizable at that level. 
The virtual, soft and collinear corrections amount to corrections in the form factors for the 
three vertices in the diagrams. The real gluon corrections may also be grouped similarly into 
three well defined channels (tt production, t decay, and t decay). This is because one may 
always identify on-shell t and t states in each diagram. Thus, one can construct a very efficient 
event generator. 

To account for gluon interference effects, one may not rely on helicity angles, because the 
amplitudes are no longer factorizable, and need to be expressed in terms of a single Lorentz 
frame. This complicates their calculation. To simplify the expressions, we have introduced 
polarization vectors for the intermediate W ± boson states defined with respect to the b (b) 
quarks (which are treated as massless). This quadruples the number of diagrams, but simplifies 
the final expressions. At tree level, we have checked the accuracy of our expressions numerically 
by using Schmidt's program. At loop level, we have performed a similar check in the limit 
Vi,2 — * 0. Before we can feed our expressions including gluon interference effects to a Monte 
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Carlo event generator, we need to integrate over the momentum of the gluon (virtual, or real, 
accordingly). In general, this can be done using the standard procedure implemented via the 
symbolic manipulation program FORM M. In our case, we need to first reduce the number 
of gluon momentum instertions to two, otherwise the final expressions turn out to be too 
lengthy to handle. This can easily be done if the amplitudes are expressed in terms of W boson 
polarization vectors, as described in the previous Section. After the amplitudes have thus been 
massaged, they can be straightforwardly reduced to scalar amplitudes via FORM || . The final 
expressions can then be fed to an event generator. 

Since gluon interference effects are relatively small, we have taken advantage of Schmidt's 
efficient code, by keeping the event generating procedure and introducing the new effects as a 
correction to the weight. In general, this correction is an extra factor to the weight, 

Wneu, = = 1 + Ax + A 2 + 0(vh) (51) 

There is an additional complication in the case of hard gluon emission, where it is no longer 
possible to define channels corresponding to the three main vertices in the diagrams (tt pro- 
duction, t decay, and t decay) by identifying on-shell intermediate t and t quark states. To 
maintain the structure of the program, we define the three channels by dividing the phase space 
of the final states into distinct regions corresponding to the three channels, as described in the 
previous Section. These definitions of the channels agree with Schmidt's in the limit 771,2 ~~ * 0. 

In our approach, we have set the mass of the b quark to zero in matrix elements. This 
was essential for the simplification of the expressions for the various amplitudes. We have also 
neglected effects of o(rjf 2 ), which include gluon interference between the b and b quarks. Such 
effects are very small and are not likely to contribute at the desired accuracy for the NLC. 
Should the need for higher accuracy arise, our method can be readily applied to accommodate 
these effects. We have also not coupled the program with a jet algorithm. This is necessary for 



a direct comparison with experimental data. Initial-state radiation effects |L3| are included in 
the program, but will not be discussed here. 

Next, we present numerical results obtained with the Monte Carlo event generator. As 
has been shown, there is no correction to the total cross-section, but differential cross-sections 
are altered by gluon interference effects. Our results include hard gluon effects and are in 
general agreement with the semi-classical analysis of soft gluon production of Khoze, et al. ||10|| . 
Specifically, we plot the distribution of gluon energy (Figs. |J H), transverse momentum (Figs. ||, 
0), and polar angle (Figs. |8], [5]). All these plots are for center-of-mass energy yfs = 400 GeV. 
We have also introduced the jet resolution parameter r\ cut = 0.05 by demanding 

(k + p b ) 2 /s > r, 2 cut , (k + Pi ) 2 /s > r] 2 cut (52) 

where k^ is the four-momentum of the gluon. We have imposed no other cut on the gluon. 

In Figs. |], |, and || we have compared distributions with and without interference for the 
three variables we are studying. In general, interference effects enter at a 10 % level, so they 
need to be included in a precision study of top quark production and decay. Moreover, in 
Figs. |5], [7|, and we show the dependency of the respective distributions on the decay width of 
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the top quark, T t . In Fig. |, we see that the decay width has a most pronounced effect in the 
range of gluon energies E g ~ T t @, and similarly for the gluon transverse momentum (Fig. |7|). 
Also, the effect of T t can be amplified in the distribution of the gluon polar angle, if one tags 
the top or bottom quarks J7|. 

VI. CONCLUSIONS 

We have performed a study of gluon interference effects in the tt production above thresh- 
old and the subsequent semi-leptonic decay at the NLC. Our study extends previous analytical 
analyses of soft gluon production, by including hard gluon effects. We have performed a numer- 
ical analysis which is complete to o(T t /m t ). We have obtained analytical expressions for the 
amplitudes, aided by the symboic manipulation routines based on FORM || . These expressions 
were then fed into a Monte Carlo event generator based on code written by Schmidt. 

Although there is no correction to the total cross-section, differential cross-sections are 
affected by gluon interference at the 10% level, in general. Therefore, such effects need to be 
included in a precision analysis of top quark production and decay at the NLC. A complete 
study should span the beam energy range including the region near threshold. This presents 
the additional complication of color Coulomb effects which require non-perturbative methods 
to be dealt with. The region near threshold has recently been studied in detail. We hope to 
extend the range of validity of our Monte Carlo event generator by including Coulomb effects 
at threshold in the near future. 

We have concentrated on the semi-leptonic decay channel, because this is where one obtains 
the cleanest signal. The branching ratios for other (hadronic) channels are higher and gluon 
interference effects are more pronounced there. We are planning to continue with a systematic 
study of the other top decay channels. Our method can also be applied to hadronic colliders, 
such as the LHC, but a similar study in that environment is far more involved, because of gluon 
interference between initial and final states. 
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FIG. 3. tt production and decay with gluon emission. 
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FIG. 5. Distribution of gluon energy for T t = 0.7,1.5,10 GeV (including interference) at cen- 
ter-of-mass energy of 400 GeV. 
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FIG. 6. Distribution of gluon transverse momentum with (dotted line) and without (solid line) 
interference at center-of-mass energy of 400 GeV. 
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FIG. 7. Distribution of gluon transverse momentum for F t = 0.7, 1.5, 10 GeV (including interfer- 
ence) at center-of-mass energy of 400 GeV. 
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